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Introduction 

When using a Groebner basis to solve the highly symmetric system of algebraic equations 
defining the cyclic p-roots, one has the feeling that much of the advantage of computerized 
symbolic algebra over hand calculation is lost through the fact that the symmetry is im- 
mediately "thrown out" by the calculations. In this paper, the problem of finding (for all 
relevant primes p) all cyclic p-roots of index 3 (as defined in Section 1) is treated with the 
symmetry preserved through the calculations. Once we had found the relevant formulas, 
using MAPLE and MATHEMATICA, the calculations could even be made by hand. On 
the other hand, with respect to a straightforward attack with Groebner basis, it is not even 
clear how this could be organized for a general p. 

In other terminologies, our results involve listings of all bi-unimodular sequences con- 
stant on the cosets of the group G? of cubic residues, or equivalently all circulant complex 
Hadamard matrices related to Go (cf. [3]). 

The corresponding problem for bi-unimodular sequences of index 2 was solved by the 
first named author in [2] and shortly after solved independently by de la Harpe and Jones 
[8] in the case p = 1 (mod 4) and by Munemasa and Watatani [11] in the case p = 3 
(mod 4) , see also [7], sect. 3. 

The organization of the paper should be clear from the section headings with the un- 
derstanding that "the main problem" refers to simple sequences of index 3 (cf. Definitions 
1.2, 1.3, and 1.4). 
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1 Notation, definitions, and problem formulation 



We begin by quoting from [2] and [3] definitions of and relations between bi-unimodular 
p-sequences and cyclic p-roots for any positive integer p. For any p-sequence x, that is any 
sequence x = (x , . . . , x p -i) of p complex numbers, define its normalized Fourier transform 
by x v = -j= Yuj=o XjU^, where uj = exp(^p). The sequence x is called unimodular if \xj\ = 1 
for j — 0, 1, . . . ,p — 1, and it is called bi-unimodular if both x and x are unimodular. 
Taking all indices modulo p, we define the periodic autocorrelation coefficients 7& by 

ik= ^2 Xjx j+k . (i.i) 

j (mod p) 

Then, by the Parseval relation and an easy calculation, 

x is unimodular (70 = p and 7^ = when k ^ (mod p)). (1.2) 

We will now express the property of bi-unimodularity with the help of a certain system 
of algebraic equations. Let z = (z , . . . , z p _i) G C p . We will call z a "cyclic p-root", if z 
satisfies the following system of p algebraic equations: 

z + zi H h 2 p -i = 0, 

Z Z 1 + ^1^2 + " " ' + Z p -iZq = 0, 

; (1.3) 

ZqZi ■ ■ ■ Z p ^2 + Z1Z2 ■ ■ ■ Zp-i + • • • + Zp^iZo ■ ■ ■ Z p ~3 = 0, 

Z Z 1 ' ' ' z p-l — 1- 

(Note that the sums are cyclic and contain just p terms and are in general not the elemen- 
tary symmetric functions.) Let now x G C p and z G C p be related by 

Zj = x j+1 /xj (1.4) 

(with x p := Xq). Clearly x is unimodular iff Xj = l/xj(\fj). In this case, (1.1) for k = 
1, 2, . . . ,p— 1 becomes the /c'th equation of (1.3). Let us call x normalized if x = 1. Then 
(1.2) can be expressed as follows: 

Proposition 1.1 A normalized x — (1, x±, X2, ■ ■ ■ , x p -\) is bi-unimodular if and only if the 
corresponding z is a unimodular cyclic p—root. 

In the rest of the paper, p will be a prime = 1 (mod 6), and we will define s := (p— l)/3. 
The multiplicative group Z* on Z p \{0} is cyclic (cf. [3]) and has a unique index-3 subgroup 
Go (the group of cubic residues modulo p). Let G\ and G2 be the other two cosets of Go 
in Z*. (The choice of the subscripts 1 and 2 will be specified later.). We will now for 
p-sequences define a property of index 3 meaning "taking few values in a way governed by 
G ": 
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Definition 1.2 We will say that x G C p is simple of index 3, if there are complex numbers, 
c ,Ci, and c 2 , such that 



Xj = c k when ^ j G G k (k — 0, 1, 2). (1.5) 

Note that we have slightly changed the notation from [2] where index 3 was called 
"pre-index 3" and where "index 3" excluded the case of index 1. i.e. Cq = C\ = c<i- 

Allowing shifts and multiplication by exponentials in a way familiar in Fourier transform 
theory, we make the following definition: 

Definition 1.3 We will say that x GC P has index 3, if for some fixed elements r ^ and 
I of Z p and some simple y of index 3 we have 

xj = u rj yj-i, (1.6) 

which amounts to 

Xj = cu rj c k when ^ j - I G G k (k = 0, 1, 2). (1.7) 
We will now define simple and general cyclic p-roots of index 3: 

Definition 1.4 By a cyclic p-root of index 3 we will mean a cyclic p-root z such that the 
corresponding x, as defined by (1.4) has index 3. We will also call a cyclic p-root z simple 
of index 3, if the corresponding x is simple of index 3. 

Note that we do not require x (and thus z) to be unimodular. 

The purpose of the present paper is to find explicitly all cyclic p-roots of index 3 (for 
every relevant prime p) using a method which utilizes the symmetries of the system. 

We will now show (following [2]), that if z is a simple cyclic p-root of index 3 and its 
corresponding x is normalized by xo = 1, then the system (1.3) reduces to a system of 
three equations for cq, c\ and C2- (To help the reader, an example is given at the end of 
the section.) Let g be a generator for Z* and let Go, G±, G2 be the cosets of Go, numbered 
in such a way that G k = {g fc+3m ; m — 0, 1, . . . , s — 1}. For every % and k — 0, 1, 2, and 
every d — 1, . . . ,p — 1, we define the transition number n ik (d) as the number of elements 
b in {1, 2, . . . , (p — 1)} for which b G G« and b + d G G k . (Subscripts are taken modulo 3. 
We do not count b = p — d). Suppose now that d G G a , i.e. that d = g a+3m for some m 
(congruences are modulo p). For each b which contributes to n ik (l), we have b = g t+3u and 
b + 1 = g k+3v for some u and v. Thus, from d(b + 1) = db + d we get 

gk+a+3(m+v) _ gi+a+3(m+u) _|_ ^ (18) 

Writing instead of n^l), we thus get 

n i+a ,k+a(d) = n ik . (1.9) 

Let us now consider a simple cyclic p-root of index 3, and let the corresponding x be 
normalized by xq = 1 and have values given by (1.5). Fix d such that d G G a , and 
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consider the individual products in the degree d equation of (1.3). These products will 
take the values (c^ + a)/(q + a) with the frequency n i+a ^ +a (d), the value c a /l once (since 
(p — 1) G Go), and the value l/c a once (since p — d G G a ). Thus (1.9) implies that all 
equations whose degrees d belong to the same coset G a , are identical, and the system (1.3) 
consists of the following 3 equations (where riik = riifc(l) are the transition numbers, and 
the c subscripts are counted modulo 3): 

Y + ~ + = 0, (a = 0,1,2). (1.10) 

° a k =o i=0 Q+a 

We will now return to the choice of the subscripts in G\ and G 2 . Without loss of 
generality, we can (and do in fact from now on) suppose that 

n 02 > n 01 . (1.11) 

In fact, we must have no2 ^ n m (see Corollary 2.3), and if n 02 < noi, we replace the 
generator g by g' :— g 2+3j , for some j such that 2 + 3j is relatively prime to p — 1. Since 
g G G\ and g' <E G 2 , this will interchange G\ and G2, and we have arrived at (1.11). 

Finally, we will give the promised example: Let p = 13, and take g = 2 or 11. Then 
G = {1, 5, 8, 12}, Gi = {2, 3, 10, 11}, G 2 = {4, 6, 7, 9}, and we will have n 00 = 0, n 01 = 
nw = n 12 = n 21 = n 22 = 1, and n 02 = n 20 = n n = 2. 



2 Number theoretic results used 

In this section we give some relations between the transition numbers defined in (1.9) 
and appearing in (1.10). These relations will lead to explicit formulas for the n^. 

The mapping b — > p — b from Z p to Z p will leave each one of the sets Gj invariant and 
thus we have 

Uij = riji, i,j = 0, 1,2. (2.1) 

2 

Moreover ^^Uij = |}((Gj \ {p — 1}), and thus (recall that we have defined s = 
3=0 

2 22 

j=0 j=0 i=0 

We will get one more linear relation between the in the following way: By (1.9), all 

noi(d) with d belonging to the same G a are equal. Thus, since |j(Go) = tl(Gi) = §(G 2 ) = s, 
p-i 2 

we get s ■ s = ^^noi(c?) = s • n_ aj i_ a , which becomes 

d=l a=0 

^01 + n 12 + n 2 o = s. (2.3) 
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With the help of (2.1), (2.2) and (2.3) we can express all our nine transition numbers 
in terms of n i and n 02 : 

n o = s - 1 - n i - n 2, 

nu = n 20 = n 02 , ^ ^ 

n 22 = n w = n i, 
n-12 = n 2 i = s - n i - n 02 . 

These relations are given in [2] and also in [5], Exercise 4.29 (d). There is, however, 
one further equation satisfied by the transition numbers. We first state this equation in 
terms of n 12 ,noi and n 02 : 

Proposition 2.1 Letp be a prime = 1 (mod 6), and letni 2 ,n i andn 02 be the transition 
numbers defined in Section 1. Then 

n m n Q2 + n 01 n 12 + n 02 n 12 = n 2 m + n\ 2 + n\ 2 - n 12 . 

We have proved this result by establishing the following explicit formulas for the con- 
volutions F * G (defined by (F * G)(a) = J2bez p F( a ~ ^)^(^)) °f certain complex- valued 
functions F and G on Z p . Let r, be the characteristic functions xg 3 of Gj (j — 0, 1,2), 
and let / = X{o}- Then, (with indices taken modulo 3): 

r i *r i = n M r + n i+ 2,i +2 ri + n i+ i )i+ iT 2 + si, 

Tj * Tj+i — ni^ + iT + nj + 2,jr 1 + ni + i t i +2 T 2 . 

Our original proof of Proposition 2.1 used these formulas and the commutativity and 
associativity of the convolution. Also, the reader of [5] is encouraged in Exercise 4.29 (e) 
to prove this proposition. But it turns out that Proposition 2.1 is just a reformulation of 
a theorem of Gauss (in Disquisitiones, Article 358), which we give in a form a little more 
precise than in [10] or [13] or [5]: 

Proposition 2.2 Letp be a prime = 1 (mod 6), and letni 2 ,noi andn 02 be the transition 
numbers defined in Section 1. Then there are integers A and B such that 

4p = A 2 + 27 B 2 . 

If we require that A = 1 (mod 3) and B > (which is always possible and which we 
always do), then A and B are unique, and we have 

A = 9n 12 — p — 1 and B = \n 02 — n i|. 

Since 4p is not a square, we must have B ^ 0, and hence we get the following corollary, 
which we needed at the end of Section 1: 

Corollary 2.3 Let p be a prime = 1 (mod 6), and let n m and n$ 2 be the transition 
numbers defined in Section 1. Then n i ^ n 02 . 
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Recall that we have in fact chosen G\ and G2 in such a way that n 2 > "01 • Since B > 0, 
we thus have 

A = 9ni2 — p — 1 and B = no2 — "01 • (2.5) 

Solving the linear system given by (2.4) and (2.5) for n ik , we have proved the following 
corollary of Proposition 2.2: 

Corollary 2.4 Let p be a prime = 1 (mod 6), let n ik be the transition numbers defined 
in Section 1, and let A and B be the numbers given in Proposition 2.2. Then 

n 12 = n 2 i = l(p + A + l), 
no2 = n 2 o = n 11 = ±(2p - A + 9B - 4), 

(2.6) 

^01 = ^10=^22 = ^(2p - A - 9B - 4), 
n o + "ii + "22 = |(p-4). 



Proof of Proposition 2.1: Starting from Proposition 2.2 and replacing A and B by 
the expressions given there and then replacing p by the expression p = 3(n i +n±2 + n 2 o) + 1 
from (2.3) we get 

= A 2 + 275 2 — 4p = -36(n in 2 + n ini 2 + n 02 "i2 - "01 _ n oi _ n 02 _ n i2 + "12) 

which completes the proof. 

Proof of Proposition 2.2: The calculations needed are given very explicitly in [13]. 
In fact the theorem of Gauss stated there in Section IV.2 is our Proposition 2.2 except 
that the statement of the theorem does not contain the value of B and for A gives the 
value M p — p — 1, where M p is the number of solutions (x, y, z) in T? v of x 3 + y 3 + z 3 = 
in the projective sense. In the proof of the theorem, the formula mB = [STT] — [STS] is 
given where m is our s, where R is our G , S and T are our G\ and G 2 (in some order), 
and where finally the symbol [XYZ] is defined for subsets X, Y, Z of Z p as the number of 
triples (x, y, z) such that x E X , y E Y , and z G Z and x + y + z = 0. In the course of the 
proof it is also shown that mM p = 9[RTS]. Thus all that remains for us to have a proof 
of Proposition 2.2 is to check that [G1G2G2] — [G1G2G1] = s{ri2 — ni) and [G0G2G1] = sn. 
We write x + y + z = 0asx + y = —z, and since G 2 = —G 2 , we have that 

[G i+ 2G 2 G k+ 2\ = ^2 n i+2,k+2{y) = sn ik , 

where we have used (1.9) with a = 2 and d — y. Thus [G1G2G2] — [GiG^Cn] = s(n2o — "22) 
and [G0G2G1] = sni2, and the result follows from (2.4), which completes the proof. 
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J_ 
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p— 4 
3 
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p— 4 
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J_ 
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3 Reduction of the main problem 

Let p be a prime of the form p = 3s + 1, s G N and let 

4p = A 2 + 27 B 2 

be the Gauss decomposition of 4p, i.e. A, £? e Z, A = 1 (mod 3) and B > (cf. Proposi- 
tion 2.2). Our main problem is to find all simple cyclic p-roots of index 3, i.e. to solve the 
set of equations (cf. 1.10 and Corollary 2.4) 

-02 fe + S)-noi + * 
--o 2 (^ + t)--oi(f + t) (3.1) 

p + A + l 2p-A + 95-4 2p- A-9B -4 , , 
"12 = g , n 02 = — , n m = — . (3.2) 

Proposition 3.1 Assume (co,ci,c 2 ) is a solution to (3.1). Then the numbers 

h, = C -^ + ^, j = 0,1, 2, (3.3) 

Cj+l C j+2 

( index counted modulo 3 ) are up to a cyclic permutation given by 

^a + T^cos^-y/), 7 = 0,1,2, (3.4) 

where 9 = |Arccos (^f) an( ^ ^he pair (£1,771) is one of the following 4 pairs: 

(3.5) 



with 



{ 




= 2 




if 


= o, 


SI 




p 2 -6p+2,4 
p 2 ~3p-A 






6^(p-4) 




2pA- 


-9p-4+3- v /p(p+4A+16) 



t(2) 

^1 2(pA+3p-l) 

(2) 3yp(p+2)-3pVp+4A+16 
''1 pA+3p-l ' 

.(3) _ -2pA-9p-4-3 v /p(p+4il+16) 

^1 ~~ 2(jM+3p-l) 

(3) 3yp(p+2)+3 P Vp+4A+16 
''1 pA+3p-l 



(3-6) 



(3.7) 



(3.8) 
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Remark 3.2 a) Let us first check that all the above formulas give well-defined real num- 
bers: Since p > 4 and \A\ < 2^/p we have 

p 2 - 3p - A > p 2 - 3p - 2VP = VPiVP - 2 )(Vp + I) 2 > 0. 

Moreover, 

p + AA + 16 > p - 8 VP + 16 = ( VP - 4) 2 > 

and since A = 1 (mod 3), we have \A + 3| > 1. Hence 

\pA + 3p-l\>\(A + 3)p| - l>p-l>0. 

b)We do not prove in this section that all four cases (3.5)-(3.8) actually occur. However 
this will follow from the proof of Theorem 4.1 in the next section. 

Proof of Proposition 3.1: To make our method of proof more transparent, we first 
consider the case p = 7. In this case A = B = 1, nyi = no2 = 1, and n m = 0. Put 

fj — Cj -\ and 



hj 



C 3+2 + Cj+l 
Cj+1 C j+2 



Then (3.1) becomes 



Consider now the matrix 



fo 




-1 


- h 


- h 


< fl 




-1 


-h 


- h 2 


h 

V 




-1 


- h 2 


- h 




2 


fo 


fi 


h 


K = 


fo 


2 


h 2 


hi 










fi 


h 2 


2 


h 




h 


hi 


h 


2 



Since 



K 



1 

Cl 

c 2 



r 1 1 1 . 

[1, -,-,-] + 

C Ci c 2 



CO 

J_ 

ci 
J_ 

C2 



[1,C , Ci,C 2 ], 



we get (considering K as an operator on column vectors) 



range(fT) = span < 



1 




1 


Co 




J_ 

CO 


Cl 




J_ 

Cl 


c 2 




J_ 

C2 



(3.9) 
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Hence rank(fT) < 2, and thus all 3 x 3 submatrices of K have determinant = 0. 
Let L = (£ij)fj =1 be the co-factor matrix of K, i.e. 



det(^), 

where is the 3x3 minor of K obtained by erasing the i'th row and the j'th column. 
Put 

Pi = hi 



P2 = il2+il3 + il4 
P3 = ^22+^33+^44 
PA = ^23+^34 + ^42- 



(3.10) 



Since £ij = for all i and j, we have in particular 



Pi=P2=P3= Pa 



0. 



This gives four equations of degree three in (/ , fi, f 2 , ho, hi, h 2 ), but taking (3.9) into 
account, we can consider Pi,P2,P3,Pa as polynomials in (ho, hi, h 2 ) only, namely 



Pi 



8 - 2(hl + h\ + hi) + 2h hih 2 



p 2 = 12 - A(h + hi + h 2 ) - 3(h 2 + hj + h 2 ) - A(h Q hi + hh 2 + h 2 h ) 
-(hi + h\ + hi) + 2(h hj + h x h\ + h 2 h\) + 3h hih 2 , 

p 3 = 12 - U(h + h + h 2 ) - 8(h 2 hl + h\) - 2(h hi + h x h 2 + h 2 h ) 
+2(h h 2 ± + h x h\ + h 2 hl) + A(h 2 G hi + h\h 2 + hlh )6h hih 2 , 

p A = 6 + 3(h + hi + h 2 ) + (h 2 + h\ + h\) + 5(h hi + h x h 2 + h 2 h ) 
—2(h h\ + h\h\ + h 2 h 2 §) — 6h hih 2 . 

Let si, s 2 , S3 denote the three elementary symmetric polynomials in h , hi, h 2 . 

Si = h + hi + h 2 
s 2 = h hi + h x h 2 + h 2 ho 
S3 = hoh x h 2 
and let a denote the antisymmetric polynomial: 

a = (ho - hi)(hi - h 2 )(h 2 - h ). 

Then, 



(3.11) 



(3.12) 



hi 
hi 



h\ + h 2 2 = s\-2s 2 



h\ 



h\ 



3s 1 s 2 + 3s 3 



h hj + hihj + h 2 h\ = |(sis 2 - 3s 3 + a) 
h\hi + h\h 2 + h 2 2 ho = Usis 2 — 3s 3 - a). 
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Hence Pi,P2,P3,P4 can be expressed as polynomials in si, s 2 , S3 and a. One gets 
pi = (8 - 2s?) + 4s 2 + 2s 3 

p 2 = (12 + 4si - 3s? - si) + (2 + 4 Sl )s 2 - 3s 3 + a 
p 3 = (12-14s 1 -8s?) + (14 + 3s 1 )s 2 -3s 3 -a 



Pi 



(6 + 3si + s?) + (3 - si)s 2 - 3s 3 



a. 



Therefore the equations p\ = p 2 = P3 = Pa = can be rewritten in the form 



8 -2s? 


4 


2 







1 







12 + 4si - 3s? - s? 


2 + 4s x 


-3 


1 




$2 







12 - 14si - 8s? 


14 + 3si 


-3 


-1 




S3 







6 + 3si + si 


3 — si 


-3 


-1 




a 








(3.13) 



A necessary condition for the existence of solutions to this system of equations is that the 
determinant of the coefficient matrix M is 0. One finds 

det(M) = 8(si - 6)(si + l)(s? + 9si + 15). 

Thus Si must be one of the 4 numbers 



6, 



1, 



,(2) 



S = 



-9 + V21 



or 



.(3) 



-9 - V21 



Let MW be the matrix obtained by substituting si = s^ in M (i — 0, 1,2,3). It is easy 
to compute the kernel for i — 0, 1, 2, 3. One finds dim(ker(MW) = 1 in all cases, and 
(for convenience writing vectors in row form) 

ker(M (0) ) = span{[l,12,8,0]} 
ker(M (1) ) = span{[l,-2,l,7]} 



ker(M (2) ) = span< 



1,-9 + 721, 79 -189 + 42^] } 



ker(M (3) ) = spanj 1, -9 - V21, Z^±iZ^5l > -189 + 42^ 
Hence there are exactly 4 solutions (si, s 2 , s 3 , a) to (3.13): 





= 6, 


s {0) - 

6 2 — 


12, 




s (o) _ 

b 3 — 


8, 




= 




), 




= -1, 


6 2 — 


-2, 




b 3 — 


1, 


a« 


= -7 




), 


(4 2) 


2 ' 


6 2 — 


-9 + 2^ 


/2T, 


s (2) - 

6 3 — 


79-17^ 
2 ' 


a< 2 ) 


= -189 


+ 42^ 


/2T ) 


(«?> 


-9-V2T 
2 ' 


b 2 — 


-9-2-\ 


/2l, 


s (2) - 

s 3 — 


79+17^21 
2 ' 


a( 3 ) 


= -189 


-42^ 


m ) 



(3.14) 
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However, there is a hidden relation between s±, s 2 , S3 and a, namely a 2 is a symmetric 
polynomial in (h ,hi,h 2 ) and can therefore be expressed in terms of Si,s 2 and s 2 . One 
finds 



2 2 
s 1 s 2 



4s?s 3 



4sn + 18sis 2 s 3 - 17 sl 



(3.15) 

It is elementary to check that this equality holds for each of the four sets (s±\ s 2 \ s% \ a^) 
found above. 

We must now in each case find h ,hi, h 2 by solving the 4 equations: 



ho + hi + h 2 = s 



CO 



h hi + hxh 2 + h 2 h = s 



hoh±h 2 = s 



(h - - h 2 )(h 2 - h ) 



(3.16) 



The solutions (h ,hi,h 2 ) to the first 3 equations in (3.16) are exactly the three roots (in 
arbitrary order) to the polynomial 



(3.17) 



Since (3.15) holds in each of the four cases 2 = 0,1, 2, 3, we have 

(h - - h 2 )(h 2 - ho) = ±a«. 

Hence the 4'th coordinate in the solution to the equations (3.13) only determines the cyclic 
order of the three numbers (ho, hi, h 2 ). For i — 0, (3.17) becomes 

h? - Qh 2 + 12h - 8 = 0. 

Hence h — hi — h 2 — 2 which corresponds to case (3.5) in Proposition 3.1. 

In the cases % = 1, 2, 3 we solve (3.17) by the classical trigonometric formula in the form 
of Lemma 3.5 below, where we use (3.30) when a < and (3.32) when a > 0. This will give 
the correct cyclic order of (ho, hi, h 2 ). Note that Lemma 3.5 can be applied because in all 3 
cases (i = 1, 2, 3) si, s 2 , S3 and a are all real (being solutions to the real linear system (3.13)) 
and thus a 2 > 0, which by (3.15) means that s\ s\ — 4sfs 3 — 4s;] + 18sis 2 s 3 — 27s 2 = a 2 > 0. 
Hence, up to cyclic permutation of (h , hi, h 2 ) we have 

hj =t,i + Vi cos (e - yj) 



where 



6 

m 
e 



\si 



= -sign(a) • \(s{ - 3s 2 )* 
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It turns out that 0® = -Arccos ( — ■= 

3 \2y/7 



in all three cases {i = 2, 3, 4), while 



3 
2 



This gives case (3.6), (3.7), and (3.8) respectively in Proposition 3.1 in the case p = 7. 

Consider now a general prime p, p = 1 (mod 3). This case is mathematically no more 
difficult than the case p = 7 but a computer algebra language as MAPLE or MATHE- 
MATICA is helpful for bookkeeping purpose. Using (3.2) and (3.1) instead of (3.9), the 
polynomials (3.11) again becomes polynomials in si,s 2 ,s 3 ,a, namely 



(3.18) 



Pi 




mn 


"7-12 


mi 3 


m i4 




1 


P2 




m 21 


m 22 


m 23 


m 24 




S2 


P3 




m 31 


m 32 


m 33 


m 34 




s 3 


Pa 




77l4i 


m 42 


m 43 


777 44 




a 



where the m^rs are the following 16 polynomials in s±: 

m u = -2s\ + 8 

77712 = 4 

77713 = 2 

77714 = 

Trial = \(A + p + l)s\ + \{2A- 7p + 20)sl + 4si + (4p 

m 22 = |(A + p + l)si + |(4p-2A-20) 

"723 = -A - 2 

777 2 4 = B 

m 31 



16) 



= lAp 2 - pA - 7p + a 2 + 2A + A)s? - £(pA + Up + 17)s 2 



m 32 
™33 

m 34 

77l4i 

m 42 

777 43 

m 44 



-|p(p-4) S i + |(2-p 2 + 8p) 
= ^(-6A - 12p - 3A 3 - 8 + 2pA)si + |(6p + pA + 14) 
= |(2A + A 2 -2p + 2) 
= -\{A + 2)B 

= \{7A-p 2 + Ap + 2A 2 +pA + 5)sl + jf(-QA + Qp + pA- 16)si 

+ |(P 2 - 4p - 12)si + |(p 2 - 8p + 16) 
= -^(9A + 3A 2 +pA + 8)si + |(6A-pA + 28) 
= |(2A+ A 2 -2p + 2) 
= ~\{A + 2)B. 
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Since pi = p 2 = P3 = Pa = 0, we must have det M = where M = (rriij)fj =1 . One finds 

detM = ^(si-6)g( Sl )r(si), 

where 

q( Sl ) = (p2 _ 3 p - A) Sl + (6A + 3p 2 - 18p) 

r(si) = (pA + 3p- 1)4 + (QpA + 27p + 12) Sl + (9pA + 54p - 3Q) . { > 

It is interesting that if det M is considered as a polynomial in the independent variables 
Si,p, A,B, forgetting the relation 4p = A 2 + 27B 2 , we will get an irreducible cubic poly- 
nomial instead of q(si)r(si). By Remark 3.2, p 2 — 3p — A ^ and pA + 3p — 1 ^ 0, so the 
equation det(M) = has exactly 4 solutions (counted with multiplicity), namely 

( 4 0) = 6 

(1) 18p-3p 2 -6A 
6 1 p 2 -3p-A 

(2) -6pA-27p-12+9y / p(p+4A+16) 
S l ~~ 2(pA+3p -l) 

(3) -6pA-27p-12-9-y/p(p+4A+16) 
S l ~~ 2(pA+3p-l) 

Let be the 4 x 4-matrix obtained by substituting si = in M. We next compute 
the kernel for in each of the four cases. Let M-^ be the 3x3 minor of obtained 
by erasing the j'th row and the fc'th column of Then 

2B 

det(M 1 ( i ) ) = —=r(A + p + l)((pA + A + Ap)sf + 3pA - QA + Up). 

— I 



(3.20) 



In particular 



det(M 1 ( ?) = -^-p(p + A + l)(A + A), 
det(M 1 ( 1 1) ) = -^p( p + A + l)(Ap-A 2 ), 



det(M 1 ( i ) ) •det(M 1 ( i ) ; 



4B 2 (p + A + 1) 2 (A + 4)(4p - A 2 ) 



9(pA + 3p-l) 

Since A=l (mod 3), we have A + 4 ^ 0. Moreover 4p - A 2 = 27 B 2 > and p + A + 1 > 
(v 7 ^- X ) 2 ^ °- Hence det(M 1 ( ?) ^ in all 4 cases. Together with det(AfW) = 0, this shows 
that for all i, s, has rank 3 and thus 

dim(Ker(M w ) = 1, i = 0, 1,2,3. 
Hence in each case (i = 0,1,2,3), s^,4 i} and 

are uniquely determined by (3.18). 
Applying Cramer's rule to the last three equations in (3.18) we get 

(i) = det Afg (i) = det ilfg (i) = det Mff 
S ' detMff' 53 detM^' ° det M$ ' 
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For % — 0, (sj°\ s 2 °\ sjj ^, a^) = (6, 12, 8, 0) as in the case p = 7 and for % — 1 we have 



'i 



18p-3p 2 -6A 
p 2 — 3p— A 

g 4p 2 A-24pA+4A 2 +p 4 -21p 3 +108p 2 -144p 
J (p 2 -3p-A) 2 

20p 2 A-96pA+8A 2 -p 4 +4p 3 -360p 2 +864p 
(p 2 -3p-A) 2 

729p(p-4) 3 B 
p 2 -3p-A) 3 • 



(3.21) 



For i — 2,3, it is more convenient to express the solutions in terms of u — y/p and 
v = \/p + 4:A + 16. We get 



and 



,(2) 
'1 

,(2) 



.(2) 



,(2) 



,(3) 



,(3) 
'2 

,(3) 



,(3) 



= -3 



m 2 +mm— 4 
u 2 +uv+2 



o (m 2 +um+6m-4)(u 2 +mm-6m-4) 
6 (u 2 +uv+2) 2 

m 4 +2m 3 m-176m 2 +m : V+40mm-32 



(3.22) 



5832 



(m 2 +mm+2) 2 
Bu 2 



-3 



2 +™+2) 3 



-uv—4 



u 2 —uv+2 



= 3 



(u 2 — uv+6u— 4) (u 2 —uv— 6m— 4) 



(m 2 -mm+2) 2 
u i -2u 3 v-176u 2 +v 2 u 2 -40uv-32 



(3.23) 



5832, 



(m 2 -mm+2) 2 
B« 2 



' (u 2 -uv+2) 3 ' 

Note that all the numbers are well-defined because by Remark 3.2, p 2 — 3p — A > 0, 
p + AA + 16 > and 

{u 2 + uv + 2){u 2 -uv + 2) = -A(pA + 3p - 1) ^ 0. 

For i = 0, we get as for p = 7 that ho — hi — hi — 2 which corresponds to (3.5) in 
Proposition 3.1. It is easy to check that the identity (3.15) is satisfied for the above sets 
(si \ s%\ Sg\ a^), so as in the case p = 7 we can determine h , hi, h 2 by Lemma 3.5 where 
we use (3.30) when a® < and (3.32), when > to obtain the correct cyclic ordering. 
Note that < 0, > and sign(a( 3 )) = sign(w 2 -uv + 2) = -sign(pA + 3p - 1). We 
obtain 



hj = + V? cos ( 39^ -^-j), j = 0, 1, 2, 



2vr 



where 0^ 



|Arccos (^p) m a ^ three cases (i 



2,3,4), while 



p 2 -6p+2A 6y^(p-4) 



3p-A 



u 2 +uv—4 



-3p-A 
12m 



1 Kf , ,.)! 3 ») = (-fc 



//- + /<r+2 ' • • 2) 



(3.24) 



uv—4 
uv+2 ' 



12m 



'_t 2 — MM+2 
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Using u = y/p, v = \/p + AA + 16, we get (3.6), (3.7) and (3.8) in Proposition 3.1. This 
completes the proof of Proposition 3.1, 

Remark 3.3 It easily follows from the proof that if c = (c ,ci,c 2 ) is a solution to the 
system (3.1) and two q are equal, then they are all equal. In fact, if e.g. c\ = c 2 , then 
with h as in (3.3) we get hi = h 2 , which leads to a = 0. But since B ^ 0, it follows from 
(3.21), (3.22), and (3.23) that a ^ in all cases except the case where all hi = 2. 

Remark 3.4 (a) At a first glance it is surprising that the angle 9 in the solution formula 
above is the same for i = 1,2,3. However, this fact has a fairly simple explanation: 
Computing the linear combination 

(p - l)pi - 2p 2 -P3- 2^4 
of the polynomials pi = Pi(si, s 2 , s 3 , a) given by (3.18) one gets 

4P ~ 7 A (2si - 9sis 2 + 27s 3 ) + ABa. 

Since pi = p 2 = p 3 = P\ = and B 2 = 4p ~ 7 A , we have the following identity 

B(2s\ - 9s lS2 + 27s 3 ) + Aa = 0. (3.25) 

But if hj — ^ + r]i cos (6 — ^-j), j = 0, 1, 2, and s\, s 2 , S3, a are defined as in (3.11) and 
(3.12) one finds 

27 

2s\ - 9sis 2 + 27s 3 = — r]l cos 3^ 

and 

a = r?i sin 30. 

4 

Hence, when r\\ 7^ 0, (3.25) is equivalent to 

3V3B cos 39 — A sin 3$ = 0. 
This has a unique solution 9 E (0, |) , namely 

^lArccotf^UiArccosf^ 
3 \3VSBJ 3 \2^p 

(b) It is interesting to compare the solutions in Proposition 3.1 with the Gaussian cubic 
sum 

U v 

It is known that (cf. [9] or Section IV. 2 of [13]) that for p prime, p = 1 (mod 3), G is a 
solution to the cubic equation 

x 3 — 3px — pA = 0. 
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This equation has the 3 solutions 

^• = 2Vpcos^-y^, j = 0,1,2 

where 9 = |Arccos (^=) as * n P ro P osrt i° n 3.1. 

It is a famous problem (the Problem of Kummer) to decide for each p which of the 3 
solutions is equal to G (cf.[9] and Section 9.12 of [10] or Section IV. 2 of [13]). 

We conclude this section by stating as a lemma the classical trigonometric solution of 
a cubic equation with 3 real roots. For completeness, we recall an elementary proof (cf e.g. 
§47 of [6]) 

Lemma 3.5 Consider the cubic equation 

h 3 - s x h 2 + s 2 h - s 3 = (3.26) 
and assume that si,s 2 ,s 3 G R and 

44 - 4s ? s 3 - 4s^ + 18sis 2 S3 - 27sg > 0. (3.27) 

Then 

sj-3s 2 > 0, (3.28) 
|2sf-9sis 2 + 27s 3 | < 2(s2_3s 2 )l. (3.29) 

Moreover (3.26) has 3 different real solutions. Listed in decreasing order ho > h± > h 2 , the 
solutions are 

hj = j + l(sl-3s 2 )^ cos (e-^j, J = 0,1, 2, (3.30) 

where 

I = lArccos ( M^£lfi±|ZfiV (3.31) 
and listed in increasing order h' Q < h[ < h 2 , the solutions are 

h'^ j + 1(4- 3s 2 )* cos ^9'-^ , j = 0,1,2 (3.32) 



where 

ff = lArccos /_M^£l£l±|Z£i 



3 3 



2( S ?-3s 2 )f 



(3.33) 
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Proof: Let h , hi, h 2 be the solutions to (3.26) and define a by (3.12). Then from (3.15) 
follows that (the discriminant) a 2 > 0. Hence ho, hi, h 2 are real and different (since if e.g. 
hi = c+id, h 2 = c— id with d ^ and ho G R we would have a 2 = —4d 2 ((ho — c) 2 + d 2 ) 2 < 0, 
whereas e.g. h = hi would imply that a = ). Substituting h = u + y in equation (3.26) 
we get 

u 3 + ru + q = (3.34) 

where r = — |(sf — 3s 2 ) and q = — ^(2sf — 9siS 2 + 27s 3 ). Applying (3.15) with Si = 0, s 2 = 
r, s 3 = — g we get a 2 = — 4r 3 — 27g 2 . Since the transformation from h to u is a translation, 
the discriminant does not change and thus (3.27) becomes -4r 3 - 27g 2 > 0. Thus r < 0, 
which is (3.28). Next we consider (3.29). Squaring this relation and introducing r and q 
we give it the form | — 27g| 2 < 4(— 3r) 3 , which we have just seen is true. 
Taking u = mz in (3.34) we get the equation 

z 3 + ^z + ^- = 0. (3.35) 

We now start from the trigonometric identity cos 39 = 4 cos 3 6* — 3COS6 1 . Writing z = 
cos 9 we give it the form 

3 1 

z 3 - -z - -cos 3^ = 0, (3.36) 

4 4 

which clearly has the solutions 

Zj = cos [6 - — ) , j = 0, 1, 2. (3.37) 



4r 

We see that equation (3.35) will be identical with (3.36) if m — d— — and 
—27q 

cos 3^ = . Returning to the variable h we see that (3.37) will lead to the solu- 

V— 27r 3 

tions (3.30) to the original equation (3.26) and that we can choose 9 as in (3.31). 
Since 9 e (0, |), we have 

-1 < cos (o - ^pj < -\ < cos (o - ^pj < \ < cos9 < 1. 

Hence ho > h\ > hi- Finally, note that with the notation from (3.33) we have 0' — | — 9 
and therefore h' Q — h 2 , hi — h[, h' 2 = h and thus h' < h[ < h' 2 . 

4 Solution of the main problem 

Theorem 4.1 The set of equations (3.1) has exactly 20 solutions in C 3 . The first two 
solutions are the \- solutions" : 

2-p± yjg - 4) 
c = ci = c 2 = . (4.1) 
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The remaining 18 solutions can be obtained from the three solutions listed below by the six 
transformations 

(co,Ci,C2) — > (Cfc, Cfc + i, Ck + 2) 



(--■-) 



lC ,Ci,C 2 J 

.Cfc Cfc+i Cfc +2 < 



where k = 0, 1, 2 and indices are computed modulo 3. Pni u = y/p, v = \Jp + AA + 16 and 
6 1 = |Arccos(2^=)- T/ze three solutions are c^) = (c$\ c^\ c^) , i = 1,2,3, where 

... 27T 27T 

C W = a « + /3« cos(# - — j) + 7 (i) sin(# - — j) (4.2) 
J 3 3 

and 



a (l) 

7 (1) 
a (2) 

7 (2) 
a (3) 

7 (3) 



2 p 2 -3p-A 1 * 2 

1 yp( P -4)(A+2) 

2 p 2 ~3p-A 



— I 



p 2 —3p~A 
3y/3 Vp=Z(p-2)B 



p 2 -3p-A 



2 p 2 -3p-A 



p 2 -3p-A 



1« 2 



>__4 _|_ i u\/A+u— v^JA— u+v 

2u 2 +uv+2 ' 2 u 2 +uv+2 

(A+2)m , j (u 2 +uv+A)^A+u-v^A-u+v 
u 2 +uv+2 ' 4 w 2 +itv+2 



+ 1 
/i 



u 2 +uv+2 

1 u 2 +itf — 4 
2u 2 -to+2 



(m 2 -m«-4)Vm+h+4\/m+i'-4 

i M 2 +MW+2 



m \Zm+i;+4Vm+i'— 4 
2 u 2 -uu+2 



(A+2)m 1 (m 2 — uv+A)^u+v+4 x /u+v-4 : 



u 2 —uv+2 
3\/3Bu 



u —uv+2 



. l (M 2 +TO-4)V4+M-tV4-M+i; 
2 -to+2 4 m 2 -mw+2 



(4.3) 



(4.4) 



(4.5) 



T/ie solutions (4.2) gruen fry (4.3) and (4.4) are unimodular while the e-solutions and the 
solution (4.2) given by (4.5) are real. Hence of the 20 solutions 12 are unimodular and 8 
are rea/. 



Remark 4.2 

(a) Of course the choice of a "canonical" solution among six possible ones is arbitrary. Our 
choice is motivated by a wish to give the asymptotic results in Section 6 a simple form. 

(b) It follows from the proof of Theorem 4.1 that the transformation 

(co, Ci, C2) — > can be obtained just by changing the sign of the second term in 

the above formulas for aW, (3^\ and 7W . 
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(c) Since u = y/p and v = \pp + AA + 16 and \A\ < we have 

|u - 4| < n < -a + 4 

which means that the numbers u, v,A can be the lengths of the three sides in a non- 
degenerate triangle. Hence the 4 square roots 

\/u + v + 4, V - " + f — 4, V4 + -a — n, \/4 — -a + n 

are well defined and strictly positive. Note also that 

v 2 — u 2 — 16 . 

A = (4.6) 

and 

1 /l 7^ V u + v + A\Ju + v - 4V '4 + u - v\/ 4 - u + v . A 

B = vi^^= ■ (47) 

The proof of Theorem 4.1 relies on Proposition 3.1 and the following 3 lemmas: 

Lemma 4.3 Let a , a±, a 2 G C and let G K. Tnen tnere are unique numbers p, a, r G C 
sncn i/iai 

aj = p + a cos ^0 - + r sin ^0 - , j = 0, 1, 2. 

Proof: By an elementary computation one finds 

det I 1 cos(0- f) sin(0- f N " 3v " 
\ 1 cos(0 - f ) sin(0 - f 

In particular the determinant is non-zero, which proves Lemma 4.3. 

Lemma 4.4 Let G R and Zei «i, 71, «2, /?2, 72 G C, ana" pwt 

ai + a 2 , Pi + Pi f a 2tt A ,71+72 . /. 2tt . 

c 3 = —^ + ^™{ d -Y 3 ) + —^ sm { e -T 3 

ai-a 2 Pi-fo f a 2tt A 71 - 72 . ( a 2tt . 

c j = — 2 — + — 2 — cos V ~ T J ' — s — sm ' ~ ~ J 



/or j = 0, 1,2. Then the following two conditions are equivalent 

(i) c c = C1C1 = c 2 c 2 = I, 

(ii) t 1 =t 2 = t 3 = 0, 
where 

tl = (al-aD + ^-PD + ^-^D-A, (4.8) 
t 2 = 2(aiPi - a 2 p 2 ) + \{p 2 x -Pi -7i 2 + 72 2 )cos30 + (A7i-/3 2 72)sin30, (4.9) 
t 3 = 2(a l7l - a 272 ) + HPl ~ Pi ~ ll + 7 2 2 ) sin 30 - (/3 l7l - ^72) cos 30. (4.10) 
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Proof: Put 



fj = °j + °j = a i + A cos (° ~ yi) + 7i sin (0 - y j 
& = Cj - Cj = a 2 + /3 2 cos (e - yj) + 72 sin - yj 



Then (i) is equivalent to 



0,1,2. 



By expressing cos 2 y?, sin 2 </?, cosy? sin ip in terms of cos2<p, sm2ip (ip — 6 — ^-j) one finds 



2 , PI + l\ 



\- 2«i/3i cos - yj) + 2«i7i sin - yj) 
cos (20 - yj) + /3 l7l sin (2^ - yj) . 



Using ^j 



~j (mod 27r) one gets 



cos 29 



-fj) 



cos 3$ cos 



+ sin 3$ sin 



Hence 



where 



sin (20 - f j) = sin 30 cos (0 - f j) - cos 30 sin (0 
/? = pi + (7i cos (0 - yj) + n sin (0 - yj) , 



(0 - 2 fj) 



2n 
3 ^ 



Pi = «? + |(/3i 2 



-/ 2 1 
7i J 



01 
Tl 



2«i/3i + 



0t± 
2 

:)2 „,2 



cos 36 1 + /3i7i sin 36 1 



Similarly 



where 



2«i7i + sin 36 - A7i cos 
#J = p 2 + ^2 cos (0 - yj) + r 2 sin (0 - yj) , 



p 2 = «i + |(/?| + tI) 

cr 2 = 2a 2/ 3 2 + 2a 2 /3 2 cos 30) + /3 2 7 2 sin 30 
r 2 = 2a 2 7 2 + 2a 2 f3 2 sin 30 - /3 272 cos 30. 

Since the coefficients in the decomposition 

2tt 



fj ~~ = G°i ~ ft) + - ^2) cos ( - yi) + ( r i ~ r 2) sin (# - yi) 

are unique by Lemma 4.3, we have fj — g"j = 4, j = 0, 1, 2, if and only if 



Pi ~ P2 
This proves Lemma 4.4. 



' j v 3 

Ox — a 2 — 0, and T\ — r 2 



(4.11) 
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Lemma 4.5 Let 9 e K and /e£ c , ci, c 2 G C\{0}. Pit£ 

f--r- + — o -r -- ft . - 52±i + £z±i , . _ C7+2 _ c^+i 
J j - h ^ > yj - h ■> <h - ^ ' j ' 

where j = 0,1,2 (counted modulo 3). 

Let moreover a u , P v , r y v ,C iU ,ri v ,( u (v = 1,2) be the coefficients in the decompositions 



fj = a 1 + p 1 cos(0-fj)+>y 1 sm(e-fj) 

9i = a 2 + /3 2 cos(0-f j)+ 72 sin(0-f 7) 

hj = £i + »7iCOs(0-fj)+Ciffln(0-f j) 

kj = & + V2COs(e-fj)+( 2 sm(e-fj) 



Then 



and 



6 = 
m = 
Ci = 


-|(Q!iA - 
-|(ai7i - 


)-l 

- a 2 {3 2 ) 

" «272) 


6 = 

V2 = 
C2 = 


f (&71 - 
3 #(72«i - 
^(a 2 A - 


A72) 
7i« 2 ) 
tti/9 2 )- 



(4.12) 
(4.13) 

(4.14) 
(4.15) 



Proof: Clearly 
Hence 



fy? - l(fj+lfj+2 ~ 9j+l9j+2) 

kj = \{fj+i9j+2 - gj+ifj +2 ). 

By expressing 008(0 — ^ E ), sin(0 — ^), cos(0 — ^f), and sin(0 — -y) as linear combinations 
of cos and sin one gets 

/i/a = (a\ - A l1 ) - aiPi cos - a l7l sin 9 + ^ ~ 7l cos 20 + /3 l7l sin 20. (4.16) 
Using now (4.11) from the proof of Lemma 4.4, we have 

/1/2 - fo = -~M + 7i) - 3ai/3i cos - 3a l7l sin 0. 
Repeating the same argument with — yj instead of 0, we have 

Wi+2 - f- = -~M + 7i) - 3«iA cos (0 - yj) - 3ai 7 i sin (0 - yj) (4.17) 
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and in the same way we have 

gj+i9j+2 - g) = + ll) - 3a 2 /3 2 cos (o - yj) - 3a 272 sin (o - yj) . (4.18) 

By the definition of fj and gj we have 



Hence, by (4.17), (4.18), and (4.19) 

— fj+ifj+2 — 9j+i9j+2 



= 4 - f (ft + 7l 2 - fS\ - 7 2 ) - 3(aiA - a 2/ 3 2 ) cos(£ - f j) 
-3(«i7i - a 272 ) sin(0 - ^ j). 

By uniqueness of this decomposition (Lemma 4.3) we can read off the coefficients £i, rji, (i 
in (4.13) namely 



6 = 




+ 7i 2 -/3 2 2 - 7 2 2 ), 


»7i = 


-|(«i/9i 


- 012P2), 


Ci = 


-|(ai7i 


- 01212)- 



However by (4.8) in Lemma 4.4, we have 

(« 2 -« 2 2 ) + i(^-/3 2 2 ) + |( 7l 2 - 72 2 )= 4. 
Hence the above formula for £1 can be changed to 

ei = f(«?-«i)-i. 

This proves (4.14). A similar but much simpler computation gives 

% = k(fj+l9j+2 - fj+29j+l) 



(/3 27 i - A72) + — (72«i - 710:2) cos - yj) 
y(a: 2 /3i - «i/3 2 ) sin (d - yj), 



which proves (4.15). 

Proof of Theorem 4.1: Assume that (co, ci, c 2 ) is a solution to the set of equations (3.1). 
By Proposition 3.1, the numbers 

ft; = 2±H + 2±I, J = 0,l,2, 

Cj+1 c J+2 
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must be of the form 

2tt 



h, = £1 + 171 cos (O - yj) , j = 0, 1, 2, (4.20) 
where (£i, 771) is one of the four pairs 

(d^f), * = 0,1,2,3, listed in (3.5)-(3.8). For i = 0, 
we have £1 = 2 and 771 = 0. Hence ho = hi = hi = 2 which implies that Co = c\ = C2, and 
in this case the only solutions to (3.1) are the 2 "e-solutions" from [2], namely 



2-p± y/p(p ~ 4) 

C = Ci = c 2 = . 

For 7 = 1, 2, 3 we can compute the numbers Cj from (£1, 771) by Lemma 4.5. Define 
f.-r- + - n-c-- ft.-^±i + ^±i fr.-^±a_^±l 

Cj Cj Cj+i Cj+2 C J+ i C,+2 

as in Lemma 4.5, and let ck^, P u , ^y u , £;,, T) u , (^ u , v = 1,2 be the coefficients in the decom- 
positions (4.12) and (4.13). Note that by Lemma 4.3 this new definition of £i and 771 is 
consistent with (4.20). Moreover d = by (4.20). 
By (3.1) 

p - 4 p + A+1 2p-A-9B -4 2p - A - 9B - 4 
ti ~ 3 9 *i is ^+1 ^ W 



Since 



ho = £1 + 771 cos 0, 

/&i = £1 + 771^ - I cos0 + C^- sin#j, 

£1 + 771 ^ - \ cos - 7— sin 9 j , 



/i 2 = 



we have 



,. , P-4 P-K\ A + 2 V3 

fo = ( —3 3-6 J - —m cos - —B m . 



Repeating the same computation with 9 replaced by 77 — %Lj, we get that the coefficients 
«i,/3i,7i in the decomposition 

fj = ai + /3i COS (9 - yj) + 71 Sin (9 - yj) 

are given by 

p— 4 p-i f 

3 3 si 

771 (4.21) 
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Provided a\ - §(£i + 1) ^ we then get from (4.14) 



« 2 
fa 

72 



ICO- 



(4.22) 



Inserting the values r]?), i = 1, 2, 3 from (3.24) in (4.21) we find that a?-f (& + 1) ^ 
in all the cases i = 1,2,3. Hence the numbers a±, 71, a 2 , (3 2 , I2 given by (4.21) and 
(4.22) are unique up to simultaneous sign change of (0:2, P2, 72)- For i — 1,2, 



a 2 = ± 



n 7 



3(6 + 1) 



is purely imaginary, and we choose the solution with Q^a^) > (i — 1,2). For 2 = 3, a 2 
is real and we choose the solution with sign^ 3 '') — — sign(w 2 — uv + 2). It is now easy to 
compute cti, 71, a 2 , (3 2 , I2 explicitly from (3.24) in the 3 cases i = 1,2,3. One finds 



,(!) 



a = 



^ = 



7i 



(i) 



p 2 -3p-A ' 
yp( P -4)(A+2) 

3V3^(p-4)g 



a' 1 ' = 



pr—3p—A 



p 2 -3p-A 



(1) - ^l(pA-2p-2A) 
~ 1 p 2 -3p-A 



(4.23) 



a? = 
P? = 



7i 



u — uv— 4 

u 2 +uv+2 ; 

2(A+2)tt 

M 2 +TO + 2' 



(2) _ 6V3Bi 



u 2 +uv+2' 



(2) 

«2 = 

tf(2) _ ; 

/Jo — 9 



72 



(2) _ 



M 2 +MW + 2 



1 (u 2 +-fiv+4)V4+M-vv / 4-u+w 

2 u 2 +uw+2 

i (tt 2 -w+4)Vtt+i'+4VM+i'-4 



u 2 +m;+2 



(4.24) 



Since 





u 2 +uv—A 


<4 S> 


u 2 —uv+2 ' 


< /??> = 


2(A+2)u 
u 2 — uv+2 ' 


/?? 


(3) 

k 7i = 


u 2 —uv+2 ' 


(3) 
72 


= Wi+9i) 


CCi + « 2 

2 


2 



My / tt+i'+4\/M+i;— 4 
u 2 — uu+2 

1 (M 2 -tw+4)Vu+i;+4\Ai+i;-4 

2 it 2 -?iu+2 

1 (M 2 +M«-4)V4+M-«V4-M+i; 

2 u 2 -uu+2 



fa 2n \ , 7i+72 . / fl 2vr 
cos —j H — sin ( 9 



(4.25) 



(4.26) 



we obtain (4.2) with a®,P® n M given by (4.3), (4.4) and (4.5). 
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We still have to check that the (cq\ c(\ c^) given by (4.2)-(4.5) actually are solutions 
to (3.1). From Lemma 4.4 and Lemma 4.5 it follows that the only thing left to check is 
that Cj ^ 0, j = 0, 1, 2 and that 

1 ai - a 2 , /3i - /3 2 ( Q 2tt .\ , 7i - 72 . ( n 2tt 
H cos [0 — '' 



+ (4,7) 



9 2 2 V 3 / 2 V 3 

which is equivalent to checking that the numbers ti,t 2 ,t 3 listed in (4.8)-(4.10) are zero. 
Using 

on A ^>-A 2 3V3B 
cos3fl = — — , sm30 = — — = — 

it is elementary to check by MAPLE or MATHEMATICA that t x = t 2 = t 3 = in each of 
the 3 cases (4.23), (4.24) and (4.25) above. It is also possible to avoid a case by case check 
by relating t±, t 2 and t 3 to the polynomials Pi,p 2 ,P3,p^ used in the proof of Proposition 
3.1 (see Remark 4.6 below). 

Finally we have to show that we have found 20 distinct solutions: Since i]^ ^ 0, 
% — 1, 2, 3, the 3 solutions given by (4.2)-(4.5) are distinct from the two e-solutions. This 
also implies that in each of the 3 cases, the 6 solutions given by 



(cj,c j+1 ,c j+2 ) j = 0,1,2 
(-,—,—, 7=0,1,2 

V Cj ' c j+1 ' c j+2 ' 



(4.28) 



are all distinct. To check that there is no overlap between these 3 groups of 6 solutions it 
is sufficient to check that the 3 numbers s±^ = 3£^ are distinct because 

7 . , . , C 2 C Ci Ci c 2 c 

si = h + hi + h 2 = 1 1 1 1 1 

Cl c 2 C C 2 C Ci 

is invariant under the 6 transformations listed in (4.28). From (3.20) 

(i) 18p - 3p 2 - QA 



s 



1 ~~ ^2 



,(2) 
'1 

.(3) 



p 2 — 3p — A 
-6pA - 27 - 12 ± 9 v /p(p + AA + 16) 
2(pA + 2p- 1) ' 



Clearly ^ sf \ since p + AA + 16 > by Remark 3.2. Moreover sf^ and sf -* are the 
two zeros of the polynomial r from (3.19): 

r(si) = (M + 3p - l)s? + (6pA + 27p + 12)s x + (9pA + 54p - 36). 

We get 

r(s^ ; ) = 81- 



S1U o ^2p-A-A)(Ap-A 2 ) 



( P 2 -3 P - Ay 
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but 2p- A- A > 2p-2 v /p-4 = 2(y/p+l)(y/p-2) > 0, and 4p- A 2 = 27B 2 > 0. Hence 
7^ s ? anc ^ 7^ S f^- Therefore we have found altogether 2 + 3 • 6 = 20 solutions. 
By (4.27), passing from to -fa in (4.2) corresponds to a change of sign of a 2 , f3 2 and 

7 2 . Hence the 12 solutions generated by (4.3), (4.4) and the transformations (4.28) are all 
unimodular while the remaining 8 solutions clearly are real. 
This completes the proof of Theorem 4.1. 

Remark 4.6 We sketch here a different proof of t\ = t 2 = t 3 = for the values of 

ai,/3i,7i,a 2 ,/52,72 listed in (4.23)-(4.25): 

By (4.21) and (4.22), oti, 0i, 71, a 2 , (3 2 , 72 can be expressed in terms of (£1, 771) and hence 
ti,h,t 3 given by (4.8)-(4.9) can be expressed in terms of £1,771, and 9. Next we observe 
that if ^ 

hj = £1 + 771 cos (0 - y/), j = 0,1,2, 



then 



si = h + /ii + h 2 = 3£i, 

3 

s 2 = 7 lQ / ll + h 1 h 2 + h 2 h = 3£ 2 - -r} 2 , 
S3 = hoh x h 2 = £1 - |£i77i + \r}\ cos 39, 



3\/3 

a = (h -h 1 )(h 1 -h 2 )(h 2 -h ) = — 77isin36>. 

Inserting this into the 4 polynomials pi = p?(si, s 2 , S3, a) from the proof of Proposition 3.1 
and comparing these new formulas for Pi,p 2 ,P3 and p 4 with the formulas found above for 
ti,t 2 ,ts one discovers after some work that 

4(ps - Pa) 



27a 2 



4(3pi + (6 - g6 - p + A)p 2 + (26 - l)p 3 + gi + 4)p 4 

27£ia| 

4(6 +p£i + A£i + p + A - 2) ((p - l)pi - 2p 2 - Ps - 2p 4 ) 



by sol 

it follows that ti = t 2 = t 3 = in all three cases. 



and since (£[ , 77^ ), i = 1, 2, 3 were found by solving the equations pi = P2 = P3 = Pa = 0, 



5 Corollaries of the main result (Leaving the simple 

case) 

In this section we will formulate and prove various consequences of the main result; in 
particular we will identify all bi- unimodular p-sequences and cyclic p-roots of index 3. We 
will give the names: 
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Definition 5.1 We denote as the first, second and third canonical solution the solutions 
c^ l \ c^ 2 \ and c^ defined in Theorem 4.1. 

We will start by presenting all bi-unimodular p-sequence of index 3 (cf. Definition 1.3). 
Recall that u = exp(^). 

Proposition 5.2 Letp be a prime = 1 (mod 6), and let x be a bi-unimodular p-sequence 
of index 3. Then there are a complex number b of modulus one and integers r and I such 
that x is given by xi = b and Xj = b ■ uj r i • Ck when ^ j — I G Gk (k — 0, 1,2), 
where c = (co, c±, C2) is one of the 12 solutions to (3.1) coming from the the first or second 
canonical solution c^ l \ c^, as described in Theorem 4.1. If p ^ 7, there are Yip 2 different 
normalized bi-unimodular p- sequences of index 3 (i.e. withxo = 1). There are 336 different 
normalized bi-unimodular 7 -sequences. Of these, 6 • 7 2 come from the second canonical 
solution, whereas only 6 ■ 7 come from the first canonical solution. The last-mentioned 
sequences can be uniquely written in the form Xj = uj m '^ +n \ where m and n G Z 7 and 

Next we formulate our result as a theorem bearing on cyclic p-roots rather than on 
bi-unimodular p-sequences: 

Proposition 5.3 Let p be a prime = 1 (mod 6), and let z = (zq, ■ ■ ■ ,z p -i) be a cyclic 
p-root of index 3. Then there are integers r and I such that z is given by Zj = u) r ■ Ck/c K 
when j + 1 — / G Gk and j — I G G K , where c = (cq, Ci, c 2 ) is one of the 20 solutions to (3.1) 
as described in Theorem 4.1. If p 7^ 7, there are 20p 2 different cyclic p-roots of index 3, 
(2p 2 of which being in fact of index 1 ). There are only 434 different cyclic 7 -roots of index 
3. Of these, 42 come from the first canonical solution. These "Gaussian" cyclic 7 -roots 
can be uniquely written in the form Zj = oj m ^ +n where m and n G Z 7 and m^O. 

Proof of Proposition 5.2 and Proposition 5.3 The first statements in these theorems 
are obvious reformulations of Theorem 4.1 in terms of the concepts introduced in Section 
1, and we leave it to the reader to check this. We will only prove the statements about the 
number of different normalized bi-unimodular sequences of index 3 (NBUS3), the number 
of different cyclic p-roots of index 3, and the explicit forms given in the first canonical case 
for p = 7. 

We start with the last topic. Since the 42 possible w-exponents in the ^-formula in 
Proposition 5.3 form the set of all differences (as functions of j) of those in the Xj-formula 
in Proposition 5.2, it suffices to consider the latter (cf. (1.4) and Proposition 1.1). We 
start by taking m = 1 and n — 0, which gives x = (1, u, w 4 , u 2 , u 2 , cj 4 , u). Since for p = 7 
we have Go = {1,6}, G\ = {3,4}, and G2 = {2,5}, this means that this particular x is 
in fact simple of index 3 with Co = u>, c\ = u 2 , and C2 = oj a (cf. Definition 1.2). We 
claim that this c = (c , ci, c 2 ) is one of the six solutions coming from in Theorem 4.1. 
To prove this, we calculate h = a + a = u 2 + u~ 2 , hi = ^ + 21 = u 3 + u;" 3 , and 

f ' U C2 ci ' 1 c c 2 ' 

^2 = Co/ci + c\ /cq = uj + uj^ 1 . Thus, using the relation 1 + uj + u 2 + u 3 + u A + u 5 + u 6 = 0, 
we get si = h + hi + h 2 — — 1, s 2 = h h\ + h x h 2 + h 2 h = -2, s 3 = h hih 2 = 1, and 
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a = (hi — h )(h 2 — hi) (ho — h 02 ) = —7. Since these values agree with those of s[' , s 2 , S3 , 
and in (3.14), our last claim is proved. 

Next we keep n = but consider a general m. But all we have used about uj in 
our calculations is that u is a primitive seventh root of unity. So is uj m . Thus, Xj = 
<jj m 'i = (u m y will also give a simple bi-unimodular 7-sequence of index 3. Of course the 
six possibilities for m correspond to the six transformations mentioned in Theorem 4.1. 
Finally, taking a general n, we see by Definition 1.3 (with / = and h = n) that all our x 
are bi-unimodular 7-sequences of index 3. Clearly they are normalized. 

It is clear that the 42 normalized bi-unimodular 7-sequences of index 3 we have found 
are different. Next we show that no other normalized bi-unimodular 7-sequence comes 
from the first canonical case. All we have to prove is that taking I ^ in Definition 1.3 
does not give anything new when y is a simple bi-uninormal sequence of index 3 given by 
Uk = oj mi . But this is trivial, since Definition 1.3 gives the unnormalized bi-uninormal 
sequence x of index 3 defined by Xj = uj h 'j +rn ^~ l ) = u ml +mjl - 2ml i which is normalized 
through division by Xq = to ml and becomes uj m ^ ~ 2ml ^ = uj^~ 2rnl ^y^ which is of the desired 
form. 

It remains to prove that the numbers of different NBUS3:s and different cyclic p-roots 
of index 3 given in our two propositions are correct, that is that no such "collapse" occurs 
except in the first canonical case for for p = 7. Recall that in the end of the proof of 
Theorem 4.1 we showed that all the 20 solutions to the main problem are different. We 
now have to extend this from the simple to the general case and we start by considering 
the e-solutions. Every corresponding NBUS3 x has the form Xj = djUJ rj with r G Z and 
d — (1, e, e, . . . , e, e) or d — (. . . , 1, 1, e, 1,1,.. .) with e = (2 — p ± \Jp(p — 4) )/2. These p 2 
NBUS3:s are clearly distinct. 

Let us when r ^ and / are in Z p and c = (c , Ci, c 2 ) G C 3 is one of the 20 solutions 
mentioned in 4.1, define x(r, I, c) as the NBUS3 x = (xo, x', . . . , x p -\) given by the formulas: 

Xj = bu r:) (c)k when ^ j — I G (5.1) 

xi = boo rl , (5.2) 

where b is determined by the normalization 

x = 1. (5.3) 

Let us consider two coinciding NBUS3:s, x(r', I', d) = x(r" , I", c") which do not satisfy 
all the three equalities r' = r",l' = = c". We denote the two b:s defined by (5.3) 
by b' and b" , respectively. We start by considering the possibility that I" = I'. Denote 
the common value by I and fix a k. From (5.1) follows that b'ufi(d)k = b"uf ■ ? '( c ")fc ^ 
j — I G Gt and thus for at least two different non-zero j, which leads to r' = r" . Then (5.2) 
gives b' = b". Now (5.1) implies that we have also d = c", which is against our hypothesis 
that at least one of r, I and c differs between the two NBUS3:s. 

Thus we have V ^ I". Let us now suppose that r' = r" (and I' ^ I"). Denote the 
common r- value by r. Choose j\ such that ji ^ V and j\ 7^ I" an d define k\ and k 2 by 

( 3 i-l')eG kl , (ji-i")eG k2 . (5.4) 
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Consider the set F := {j - I"; (j - 1') G G kl } n G k , 2 . Taking d = I" - V in (1.9), we see that 
if d G G a , then the cardinality of F is a transition number: = n kl - a;k2 - a . By (2.4), 

all transition numbers are < s — 1, and since t}(Gfci) = s, there is at least one j2 and one 
^3 7^ ^2 such that 

(j2-l')eG kl and (j 2 -(")GG fe (5.5) 
Now from (5.4) and (5.5) follows that (5.1) with j = ji and with j = j 2 gives 

b'u^c' kl =b"u;^4 2 , 

This leads to d k2 = d' k , A Then it follows from from Remark 3.3 that c" is an e-solution. Since 
d and c" play the same part in our situation, the same must be true for d . But we know 
already that there is no internal collapse among the NBUS3:s coming from e-solutions, so 
the case r' = r" also leads to a contradiction. Now we know that r' ^ r" and I' ^ I" . From 
(5.2) and (5.1) with j = I' we get 

Xv = b'uo r ' 1 ' = b"u r " l 'dl, (5.6) 

where k is determined by (I" — V) G G k - Since G k has at least two elements we can choose 
j ^ I" with (j - V) G G k . For this j we get from (5.1) 

Xj = b'uj r ^d k = b"uj r J 'cfc. (5.7) 

From (5.6) and (5.7) we get by division 

d k = J r "- r 'X j - l '\ 

Since the exponent of uj is not zero (modulo p) , we have found a d i which is a primitive 
p'th root of unity. But we have also proved that we must have p = 7. For if p > 13, 
there are more than two elements in G k , and we can make two different choices of j, giving 
conflicting values to d k . To sum up, we know that to have collapse we must have p = 7, 
and some d k must be a seventh root of unity. Again our symmetry argument says that 
also some d' k must be a seventh root of unity. The third canonical case is not of interest, 
since the absolute values are not one. We can also easily exclude the second canonic case 
e.g. with the following numerical argument: The imaginary part of the seventh power of 
the six values of the components of are approximately ±0.92, ±0.94, and ±0.41 rather 
than 0. So the collapse is an internal affair within the first canonical case, which we have 
already studied. This completes the proof of the two propositions. 



6 Numerical and asymptotic results 

In this section we will study the behavior for large p of the solutions d l \i = 1, 2, 3 defined in 
Theorem 4.1. We will give numerical data leading to educated guesses about this behavior 
(see Remark 6.3 and we will prove quantitative forms of these guesses. 
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In Table 6.1 below we list the first few primes = 1 (mod 6) and corresponding numer- 
ical values of A, B, 6, c<p, and c£\ In Table 6.2, we give the corresponding information 
for . We will also include an indication of the shape of the triangle formed by the three 



complex numbers c 



(0 JO JO 

5 C l 5 C 2 



, (i — 1,2), reasoning as follows: 



In the corresponding situation for simple bi-unimodular sequences of index two (cf. [2]) 
we have two complex numbers c and C\ on the unit circle, and with increasing p their 
sum tends to zero. A natural guess in our situation might therefore be that he sum of 
the three numbers tends to zero or, equivalently, that the triangle becomes more and more 
equilateral when p grows. We prefer the latter description. To be able to give quantitative 
results we will revive the old noun scalenity, (cf. [1]) and give it a precise meaning: 

Definition 6.1 In the complex plane, let b = (b ,bi,b 2 ) be a triple of points on a circle C 
with center w. Let (pi = arg(6j — w). Let the scalenity of b be 



seal (6) 

(indices counted modulo 3). 



max 

3 



- + COs(0 j+2 



Remark 6.2 Since ^ = — cos27r/3, the triangle with vertices b will be equilateral iff its 
scalenity is zero. Let us now consider the definition of hj (in Proposition 3.1). If we 
take b = with % = 1 or 2, we have all \bj\ = 1 and thus w — 0. Hence scal(cW) = 
| maxj |1 + hj\, where hj is given by (3.4) with c replaced by b. 



Table 6.1 (First canonical case) 


P 


A 


B 


6 


Ji) 
c o 


Ji) 
c i 


Ji) 

c 2 


scal^ 1 )) 


7 


1 


1 


0.4602 


-0.9010 - 0.4339 i 


0.6235 + 0.7818 i 


-0.2225 + 0.9749 i 


1.1235 


13 


-5 


1 


0.7790 


-0.4822 - 0.8761 i 


0.3953 + 0.9185 i 


-0.8132 + 0.5820 i 


0.7132 


19 


7 


1 


0.2129 


-0.9528 - 0.3037 i 


0.9838 - 0.1791 i 


0.3780 + 0.9258 i 


0.7061 


31 


4 


2 


0.4011 


-0.8023 - 0.5969 i 


0.9923 + 0.1235 i 


-0.0963 + 0.9954 i 


0.5274 


37 


-11 


1 


0.9001 


-0.0604 - 0.9982 i 


0.4630 + 0.8863 i 


-0.9452 + 0.3265 i 


0.4127 


43 


-8 


2 


0.7423 


-0.3124 - 0.9499 i 


0.7272 + 0.6864 i 


-0.7742 + 0.6330 i 


0.3792 


61 


1 


3 


0.5022 


-0.6466 - 0.7628 i 


0.9759 + 0.2181 i 


-0.3560 + 0.9345 i 


0.3564 


67 


-5 


3 


0.6271 


-0.4569 - 0.8895 i 


0.8964 + 0.4433 i 


-0.5999 + 0.8001 i 


0.3170 


73 


7 


3 


0.3829 


-0.7843 - 0.6204 i 


0.9988 - 0.0481 i 


-0.1114 + 0.9938 i 


0.3409 


79 


-17 


1 


0.9483 


0.1286 - 0.9917 i 


0.4824 + 0.8759 i 


-0.9765 + 0.2154 i 


0.3066 


97 


19 


1 


0.0890 


-0.9875 - 0.1576 i 


0.7708 - 0.6371 i 


0.4823 + 0.8760 i 


0.3137 


103 


13 


3 


0.2919 


-0.8668 - 0.4986 i 


0.9629 - 0.2697 i 


0.0653 + 0.9979 i 


0.2937 


109 


-2 


4 


0.5556 


-0.5387 - 0.8425 i 


0.9666 + 0.2561 i 


-0.4841 + 0.8750 i 


0.2562 


127 


-20 


2 


0.8875 


0.0580 - 0.9983 i 


0.6211 + 0.7837 i 


-0.9411 + 0.3382 i 


0.2464 


139 


-23 


1 


0.9731 


0.2257 - 0.9742 i 


0.4899 + 0.8718 i 


-0.9873 + 0.1588 i 


0.2387 


151 


19 


3 


0.2290 


-0.9138 - 0.4062 i 


0.9066 - 0.4219 i 


0.1770 + 0.9842 i 


0.2452 


157 


-14 


4 


0.7212 


-0.2380 - 0.9713 i 


0.8495 + 0.5276 i 


-0.7655 + 0.6434 i 


0.2146 


163 


25 


1 


0.0683 


-0.9922 - 0.1248 i 


0.7153 - 0.6988 i 


0.4898 + 0.8718 i 


0.2411 


181 


7 


5 


0.4359 


-0.6932 - 0.7207 i 


0.9996 - 0.0270 i 


-0.2649 + 0.9643 i 


0.2092 



30 



Table 6.2 (Second canonical case) 


P 


A 


B 


6 


i?\ 
c a 




r ( ' 
c 2 


scal(c( 2 )) 


7 


1 


1 


0.4602 


0.8173 + 0.5762 i 


-0.3890 + 0.9212 i 


0.2804 - 0.9599 i 


0.7129 


13 


-5 


1 


0.7790 


0.2469 + 0.9690 i 


-0.7728 + 0.6346 i 


0.6315 - 0.7754 i 


0.9242 


19 


7 


1 


0.2129 


0.9520 + 0.3061 i 


-0.4274 + 0.9041 i 


0.0041 - 1.0000 i 


0.4058 


31 


4 


2 


0.4011 


0.8025 + 0.5967 i 


-0.6855 + 0.7281 i 


0.2171 - 0.9761 i 


0.3844 


37 


-11 


1 


0.9001 


-0.2907 + 0.9568 i 


-0.9939 - 0.1106 i 


0.8980 - 0.4399 i 


0.6830 


43 


-8 


2 


0.7423 


0.1847 + 0.9828 i 


-0.9804 + 0.1971 i 


0.7022 - 0.7120 i 


0.5126 


61 


1 


3 


0.5022 


0.6436 + 0.7653 i 


-0.8671 + 0.4982 i 


0.3730 - 0.9278 i 


0.3232 


67 


-5 


3 


0.6271 


0.4178 + 0.9085 i 


-0.9589 + 0.2837 i 


0.5632 - 0.8263 i 


0.3571 


73 


7 


3 


0.3829 


0.7932 + 0.6089 i 


-0.7822 + 0.6231 i 


0.1987 - 0.9801 i 


0.2661 


79 


-17 


1 


0.9483 


-0.4295 + 0.9031 i 


-0.8761 - 0.4821 i 


0.9657 - 0.2596 i 


0.4408 


97 


19 


1 


0.0890 


0.9887 + 0.1498 i 


-0.4791 + 0.8778 i 


-0.2383 - 0.9712 i 


0.2383 


103 


13 


3 


0.2919 


0.8757 + 0.4828 i 


-0.7239 + 0.6899 i 


0.0609 - 0.9981 i 


0.2327 


109 


-2 


4 


0.5556 


0.5288 + 0.8487 i 


-0.9508 + 0.3097 i 


0.4751 - 0.8799 i 


0.2600 


127 


-20 


2 


0.8875 


-0.2473 + 0.9689 i 


-0.9083 - 0.4183 i 


0.9253 - 0.3794 i 


0.3193 


139 


-23 


1 


0.9731 


-0.4663 + 0.8846 i 


-0.7821 - 0.6231 i 


0.9836 - 0.1804 i 


0.3135 


151 


19 


3 


0.2290 


0.9203 + 0.3913 i 


-0.6848 + 0.7287 i 


-0.0415 - 0.9991 i 


0.1997 


157 


-14 


4 


0.7212 


0.1708 + 0.9853 i 


-0.9969 - 0.0789 i 


0.7371 - 0.6757 i 


0.2520 


163 


25 


1 


0.0683 


0.9929 + 0.1190 i 


-0.4865 + 0.8737 i 


-0.2912 - 0.9567 i 


0.1941 


181 


7 


5 


0.4359 


0.7010 + 0.7132 i 


-0.9002 + 0.4355 i 


0.2959 - 0.9552 i 


0.1824 



We present the corresponding values for in Table 6.3. Since these values are real, 
we will save some space and we use this for giving the information also in another form, 

c (3) 

namely which should shed some light on the surprising behaviour of the components. 



Table 6.3 (Third canonical case) 


P 


A 


B 


6 


3 
c 


3 

c\ 




c^/VP 


c^/Vp 




7 


1 


1 


0.4602 


-1.2221 


9.4127 


2.7389 


-0.4619 


3.5577 


1.0352 


13 


-5 


1 


0.7790 


-1.4201 


-14.6415 


2.1601 


-0.3939 


-4.0608 


0.5991 


19 


7 


1 


0.2129 


-2.2521 


8.4655 


4.8488 


-0.5167 


1.9421 


1.112 


31 


4 


2 


0.4011 


-2.8168 


17.2938 


4.6888 


-0.5059 


3.1061 


0.8421 


37 


-11 


1 


0.9001 


-3.0328 


-7.1015 


2.8445 


-0.4986 


-1.1675 


0.4676 


43 


-8 


2 


0.7423 


-3.2558 


-20.3776 


3.6527 


-0.4965 


-3.1075 


0.557 


61 


1 


3 


0.5022 


-4.0014 


50.9574 


5.4586 


-0.5123 


6.5244 


0.6989 


67 


-5 


3 


0.6271 


-4.2289 


-95.9688 


5.0005 


-0.5166 


-11.7245 


0.6109 


73 


7 


3 


0.3829 


-4.4100 


25.6091 


6.6407 


-0.5162 


2.9973 


0.7772 


79 


-17 


1 


0.9483 


-5.6126 


-8.9422 


4.1623 


-0.6315 


-1.0061 


0.4683 


97 


19 


1 


0.0890 


-5.0982 


13.5365 


10.4124 


-0.5176 


1.3744 


1.0572 


103 


13 


3 


0.2919 


-5.2556 


21.9500 


8.4142 


-0.5179 


2.1628 


0.8291 


109 


-2 


4 


0.5556 


-5.5068 


337.8101 


6.6180 


-0.5275 


32.3563 


0.6339 


127 


-20 


2 


0.8875 


-7.1148 


-14.4873 


5.5161 


-0.6313 


-1.2855 


0.4895 


139 


-23 


1 


0.9731 


-8.4417 


-11.5986 


5.6060 


-0.7160 


-0.9838 


0.4755 


151 


19 


3 


0.2209 


-6.3543 


22.1314 


10.5843 


-0.5171 


1.8010 


0.8613 


157 


-14 


4 


0.7212 


-7.1235 


-36.5459 


6.8056 


-0.5685 


-2.9167 


0.5431 


163 


25 


1 


0.0683 


-6.5753 


16.4057 


13.3294 


-0.5150 


1.2850 


1.0440 


181 


7 


5 


0.4359 


-7.0841 


58.2887 


9.2448 


-0.5266 


4.3326 


0.6872 



Our observations are summarized in the following remark: 
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Remark 6.3 Our numerical observations and our results are of five kinds: 

(1) For each large p, the first and second canonical solutions are approximately symmetric 
to each other w.r.t. the origin. 

(2) Even though two large primes may be close to each other without their canonical 
solutions being close, large primes with approximately the same 9 will have approximately 
the same first canonical solutions and approximately the same second canonical solutions 
(even if the primes are not close to each other). 

(3) For large p, the first and second canonical solution each forms an approximately equi- 
lateral triangle. 

(4) For large p, the approximate positions of the nearly equilateral triangles are simple 
functions of 9. 

(5) If p is large, then all components of \c^\ are large. If in addition \A\ is small, that is 
if 9 is close to 7r/6, then |c^| is very large. 

To make it easier to guess quantitative results (making "approximately" more precise in 
Remark 6.3) we present a few more numerical results in Table 6.4. 



Table 6.4 (Large primes, close in size vs. close in 9- value) 



p 


A 


B 


e 


arg(4 ') 


29- tc 


arg( c o ) 


26* 


scal(cW) 


1003273 


973 


337 


0.354542 


-2.43320 


-2.43251 


0.70803 


0.709084 


0.002810 


1003279 


1993 


39 


0.033775 


-3.07411 


-3.07404 


0.06742 


0.067555 


0.002995 


100205473 


9733 


3367 


0.354372 


-2.43292 


-2.43285 


0.70864 


0.708744 


0.000281 



From Table 6.4 it seems that "approximately" means agreement in approximately | 
decimals. Thus quantitative results in terms of 0(^=) might seem plausible. In our quan- 
titative results we will use the maximum norm to measure distances in C 3 . We will also 
need a name for the equilateral "limit" triangle hinted atin Remark 6.3 (4), hopefully vis- 
ible in Tables 6.2 and 6.3, and present in columns 5 and 8 of Table 6.4. Thus we make the 
following two definitions: 

Definition 6.4 Let a = (ao, 01,02) £ C 3 7 then we define \\a\\ = max(|ao|, |ai|, 1 0,2 1 ) - 

Definition 6.5 Let p be a prime = 1 (mod 6) and let 9 = |Arccos (^p)? where Ap = 

A 2 + 27 B 2 and A = 1 (mod 3). We denote by d = d{p) = (ao,^,^) the (equilateral) 
triangle for which 

d J -=exp(2i(0-^)), J =0,1, 2. 

We will now state four quantitative results for the first and second canonical cases, 
where Proposition Q.j for j = 6, . . . , 9 is of the kind (j — 5) listed in Remark 6.3. (The 
discusion of kind (5) starts after Corollary 6.11 below.) 

Proposition 6.6 Let p be a prime = 1 (mod 6), and let and be the corresponding 
first and second canonical solution. Then 

|| c (l) +c (2)|| < 

5y/p 
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Proposition 6.7 Let p' and p" be primes = 1 (mod 6), let 9' and 9" be their respective 
9 -values and let c' and c" be their respective first canonical solutions. Then 



\c'-c"\\<2\e'-e"\ + ^=+ 1 



a/p 7 v^ 77 

The same result, with the constants 3 replaced by 43- ; holds if d and c" are the respective 
second canonical solutions. 

Proposition 6.8 Let p be a prime = 1 (mod 6), and let d 1 ^ and c^ be the corresponding 
first and second canonical solution. Then 



7 

scal(c (1) ) < — — and scal(c (2) ) < 



21 



2^/p 5^/p 



Proposition 6.9 Let p be a prime = 1 (mod 6), let d 1 ^ and d 2 ^ be the corresponding 
first and second canonical solution, and let d be as in Definition 6.5. Then 

|| c (i) + d \\ < 4= and \\c^ - d\\ < 

VP 5 Vp 



Remark 6.10 The constants in these propositions are not best possible but are chosen as 
compromises to make the proofs less cumbersome. Even if we restrict our claims to hold 
only for p > M for some large M, the constants cannot always be significantly improved. 
For instance, for p = 10 10 + 279 we have ||c^ — d\\ ~ For a kind of "best possible", 

result, see Remark 6.13. 

Since to each number 9 (in the interval [0, 7r/3]) there corresponds at most one p, it 
does not make sense to consider a sequence of p:s with a common 9. However, Proposition 
6.9 obviously has the following corollary, where we have used the notation 9{jp), c^Qo) and 
d 1 ^ (p) for the values of 9 and the first and second canonical solutions corresponding to p: 

Corollary 6.11 Let 9 be a real number in the interval [0, 7r/3]. Denote by d = (d , di, d 2 ) 

( 2 J 7T \ 

the (equilateral) triangle for which dj = exp[2i(9 — )J, j = 0,1,2. Let {p n }T be a 

sequence of primes = 1 (mod 6) going to infinity in such a way that lim n ^ 00 9(p n ) = 9 . 
Then lim^oo c (1) (p n ) = -d and lim^oo c (2) (p n ) = d. 

Before proving our four Propositions we will comment item (5) of Remark 6.3. In Table 
6.5 we present some more numerical values with focus on ^-values close to 0, 7r/6 and n/3. 
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Table 6.5 (Third canonical case for large primes) 


P 


A 


B 


6 








67 521 601 729 


-2 


100 016 


0.523600 


-0.577349 


779 550.5 


0.577353 


67 544 557 351 


1 


100 033 


0.523598 


-0.577347 


194 920.0 


0.577353 


250 004 500 027 


1 000 009 


1 


0.000002 


-0.500000 


1.000007 


1.000001 


250 018 500 349 


-1 000 037 


1 


1.047196 


-0.999995 


-0.999997 


0.500000 



These and other numerical results make it plausible that "large" in item (5) of Remark 
6.3, may be specified to mean "not much smaller than \yfp~ ", but it seems difficult to find 
^-independent estimates of "convergence rate" for the third canonical case. We are now 
ready to state a proposition: 

Proposition 6.12 // {p n }T ^ s an V sequence of primes = 1 (mod 6) going to infinity, 
then (with obviuos notation) for i — 0,1, and 2, 

liminf |c '- 3)(Pw)l >0,5. (6.1) 



We remark that this proposition implies that for every normalized x — (1, x\ . . . , x p -\) G 
W of index 3 coming from the third canonical case for a large p, either all \xj\, j ^ 0, are 
large or they are all small (leaving the canonical case via the transformations mentioned 
in Theorem 4.1 and leaving the simple case via Definition 1.3). 

We will now prove our five propositions. 

Proof of Propositions 6.6 and 6.6 Proposition 6.6 follows from Proposition 6.9 via a 
straightforward application of the triangle inequality. 

Similarly, Proposition 6.7 follows from Proposition 6.9 via the triangle inequality and 
the inequality | exp(2i0') — exp(2i0")| < 2\<f>' — <j>"\. 

Proof of Proposition 6.8. From (3.4) and (3.6) and Remark 6.2 we get 

(6.2) 

Since the right-hand side of of (6.2) is a decreasing function of A and A > —2^/p, we get 

^ SCal(C }< 2(p 2 -3p + 2Vp) - 3+ 2 (V p-l) 2 - (0) 

The last member of (6.3) as a function of p is decreasing for p > 9 and takes values < 3.4 
for p = 7 and 13. This completes the proof of the first part of the proposition. 
For the second part, we will use (3.7) and the identity 

fpy 



p + AA + 16 + p + 2 j (yp\Jp + AA + 16 - p - 2 j = A(Ap + 3p - 1) 
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to give the counterpart of (6.2) the form 

3(2p+v^) 



Vpscal ( c (2))<I v / p -(ji + £( 1 )| + |^)|) 



p + 2 + ^fp^p + + 16 



We can again take A = —l^fp. The resulting expression is easily seen to be < 4 for p > 100 
and for the remaining p we enter the true value of A (given in Table 6.1) to get a maximum 
w 4.1966 for p = 37. This completes the proof. 

Remark 6.13 From (6.3) we easily get the following result: For each e with < e < 1 we 
have 

3 + e ., / 3\ 2 

seal' — ■ ■ • • 



U (c »>)<i±i. i/p> ( 2 + _L), 



which could be contrasted with the fact that for p — 10 002 900 217 we have ^fp scal(cW) ps 
3.000015. 

In the proof of Proposition 6.9 we will work with a, j3 and 7 as given in Theorem 4.1 
and p, a, and r as given in Lemma 4.3. We will use the following lemma: 

Lemma 6.14 Let b' = (&' , 6^, 6' 2 ) G C 3 and b" = (&{,', 6£) e C 3 be given by 
b>. = p' + a' cos (0 - y ^ + r' sin (o - y ^ , J = 0, 1, 2, 

6f = p" + <r" cos (9 - yj) + r" sin (9 - yj) , j = 0, 1, 2, (6.4) 
w/iere ^Glcnd p', cr', r', p", cr", r" G C, Then 

\\b' - b"\\ < \p' - p"\ + ^la' - a"\ 2 + \t' -t"\ 2 . 



The proof of Lemma 6.14 is a straightforward application of the triangle inequality, the 
Cauchy inequality, and the identity cos 2 + sin 2 = 1. 

Proof of Proposition 6.9 In Lemma 6.14 we take b' = and b" = —d (cf. Definitions 
5.1 and 6.5). Then pf = a (1 \a' = [3 (1 \t' = 7W as given in (4.3), whereas p" = 0,a" = 
— cos 3$ — 2 sin 36*, and r" = — sin 3$ + 2 cos 36*, as is easily checked by introducing these 

values in (6.4) and applying the addition theorems for sine and cosine. Since cos 3^ 



and sin 3$ 



2VP 



, Lemma 6.14 shows that for the proof of the first half of Proposition 



6.9 it only remains to check that with a^\P^\ and 7^ as in (4.3) we have 
v^l« (1) 



+ 



A + i3y/3B 



2VP 



+ p 



7 (i) + 



3^35 - iA 



2VP 



< 3. 



(6.5) 



35 



Introducing the values of (3^\ and 7^ and replacing 35 a/3 by ^/4p — A 2 we can after 
some calculation treat the first term of the left member of (6.5) as follows 



Vp\ q 



p 2 — Ap 



< 



p 2 + 2p y /p 



P(2 + y/P) 



p 2 — 3p — A y p 2 — 3p + 2^/p y (p — 3) y/p + 2 ' 



(6.6) 



where the estimate comes from the facts that the second term of (6.6) is a decreasing 
function of A and that A > —2^/p. Let us denote by Q the expression under the big root 
sign in (6.5). Since the last member of (6.6) is a decreasing function of p with a value < 1.5 
for p = 31, we can prove (6.5) for p > 31 by checking that 



Q< (3-1.5) 2 = 2.25forp>31 



(6.7) 



Treating Q in the same way as we did with first term of the left member of (6.5) we find 



Q = 



2p 3 -{A + 6)p 2 + 2Ap — (2p — A — A) x Jp i -Ap 3 



p 2 — 3p 

Using a Taylor formula with rest term we have 

4 



(6.8) 



vV - 4p 3 = p 2 (l 



P 



p 2 -2p-2 + R 3 , 



where — ^ < R3 < (since p > 31). Introducing this in (6.8) we get 

_ 2p 2 - Ap - 2A - 8 - (2p - A - A)R 3 2(p 3 - 2p 2 + Ap - 2p - 3A - 12) 



p 2 — 3p — A 



p(p 2 — 3p — A) 



(6.9) 



Since the last member of (6.9) is an increasing function of A we can estimate it with its 
value for A = 2 v /p, which is a decreasing function of p and thus not larger than its value 
for p = 31, which turns out to be ~ 2.07 in agreement with (6.7). Finally, we check 
numerically the value of ^/p\\c {1) + d\\ for p = 7, 13, and 19. We find 2.59, 2.31, and 1.91, 
which are all < 3. This completes the proof of the first half of the proposition. 

For the second part of the proof we proceed in the same way but let MATHEMATICA 
help us to get a good start, namely by telling us that defining m(p) = y/p\\c^ — <i|| we 
have m(p) < m(43) < 4.1 if p < 10000. We get after some calculation 



VP\ a 



(2)1 



2p 



< 



2p 



2+P+ VpVp + AA+16 y 2 + p + ^^-8^+16 
if p > 10000. Thus to complete the proof is enough to prove that 



< 



100 
"99" 



P 



(3^ 



A + i3^/3B 



p 



7 



(2) _ 



3^3B - iA 



< 10.17 



(6.10) 
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if p > 10000, e.g. by proving that the first term of (6.10) is < 1.07 and the second term 
is < 9.1.. This can be done as in the proof of the first part, using (4.4). Just as we 
have studied functions of A restricted to the interval \A\ < we will now with the 

help of (4.6) and (4.7) write the left member of (6.10) as a function of u and v, where 
\u — 4| < v < u + 4. Again a certain square root can be estimated with a Taylor formula. 
We leave the details to the reader. 

Proof of Proposition 6.12 Inspired by the first two rows of Table 6.5 we expect infinities 
near 9 = ir/6, and thus, to avoid zeros in the denominator, we "turn everything upside 
down". Thus we want to prove that 

lim sup t|y~ — — 2- 

n-oc \ C \ >(p n )\ 

Suppose this is not true. Then (by taking subsewuences if needed) we can find a sequence 
{pn}T of primes = 1 (mod 6) going to infinity, such that 

lim rof" ^ = l j> ( 6 - n ) 

where these limits exist (finite or +oo) and \IA > 2 for at least one j (0,1, or 2). Since 
the interval [0, 7r/3] is compact, we can by again taking a subsequence (keeping the nota- 
tion {j9 n }i°) arrange that 9 — h m »woo 9(p n ) exists. Starting from (4.5) we replace A by 
2y/p cos 3$ and B by 2 v /p sin 3$ / '\/27 '. Introducing the resulting expressions for o/ 3 \ (3^\ 
and 7® in (4.2), we get y/p/c^ as a function of p and 9, which we denote by q(p, 9). We 
now fix 9 = 9 and study q(p, 9 ) as a function of p when p — > oo. Estimating various 
square roots with a Taylor formula, we get after a considerable amount of calculation: 

lim q(p, 9 )= ( - 2 cos 9 , -2 sin(# - 7r/6) , 2 sin(^ + ir/6)) . 

p— >oo V / 

A simple continuity argument (w.r.t. 9{p n and 9q) shows that with lj from (6.11) we have 
Z = -cos# , h = -2sin(0 o - 7r/6) , h = 2 sin(# + tt/6). (6.12) 

This is a contradiction, since we have supposed that \lj\ > 2 for at least one j. We 
have thus completed the proof and also substantiated the "very large" part of item (5) of 
Remark 6.3 (take l\ from (6.12) and consider for 9 close to 7r/6). 
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